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Abstract 

The quark gluon plasma which has been observed at RHIC is a strongly interacting system and 
has been called sQGP. This is a system at high temperatures and almost zero baryon chemical 
potential. A similar system with high chemical potential and almost zero temperature may exist 
in the core of compact stars. Most likely it is also a strongly interacting system. The strong 
interactions may be partly due to non-perturbative effects, which survive after the deconfinement 
transition and which can be related with the non-vanishing gluon condensates in the sQGP. In this 
work, starting from the QCD Lagrangian we perform a gluon field decomposition in low ("soft") 
and high ("hard") momentum components, we make a mean field approximation for the hard 
gluons and take the matrix elements of the soft gluon fields in the plasma. The latter are related 
to the condensates of dimension two and four. With these approximations we derive an analytical 
expression for the equation of state, which is compared to the MIT bag model one. The effect of 
the condensates is to soften the equation of state whereas the hard gluons significantly increase 
the energy density and the pressure. 
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I. INTRODUCTION 



One of the most interesting results of the RHIC program is the discovery of an extremely 
hot and dense state of matter made of quarks and gluons in a deconfined phase and which 
behaves like an ideal fluid While the production of such a plasma of quarks and gluons 
had been predicted, it was a surprise to find that this system is strongly interacting and 
very different from the originally expected gas of almost non-interacting quarks and glu- 
ons, described by perturbative QCD. This state has been called strongly interacting quark 
gluon plasma (sQGP) and there are many approaches to study its properties. The most 
fundamental approach is provided by lattice QCD Since lattice QCD has still some lim- 
itations, such as the difficulty in dealing with systems with large baryon chemical potential, 



there are several models (see, for example. 



-|5|) which incorporate the essential features of 

i the 



the full theory and which can be employed to study the sQGP. In some of them 
sQGP is treated as a gas of quasi-particles, in which the quarks and gluons have an effective 
mass. In some works, such as in 13, 6|, the sQGP was treated with semi-classical methods. 

n 

In p[ the color charges were assumed to be large and classical obeying Wong equations of 
motion. In this approach the quantum effects in the QGP are basically reduced to generate 
thermal-like masses and cause the effective coupling to run to larger values at smaller values 
of the temperature. 

The medium created in heavy ion collisions has high temperature and zero baryon chem- 
ical potential. On the other corner of the phase space, we find the QGP at zero temperature 
and high baryon number. Presumably, this kind of system exists in the core of dense stars. 
This cold QGP has a richer phase structure and at high enough chemical potential we may 
have a color superconducting phase. Because of the limitations of lattice calculations in 
this domain and also because of the lack of experimental information, the cold QGP is less 
known than the hot QGP. Nevertheless it is quite possible that it shares some features with 
the hot plasma, being also a strongly interacting and semi-classical system. 

In this work we shall study the non-perturbative effects in the cold QGP generated by the 
residual dimension two and dimension four condensates, using a mean field approximation. 

In the vacuum, non-perturbative effects have been successfully understood in terms of the 
QCD condensates, i.e., vacuum expectation values of quark and gluon "soft" (low momen- 
tum) fields. The best known are the dimension three quark condensate and the dimension 
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four gluon condensate These condensates can, in principle, be computed in lattice QCD 
or with the help of models. In practice, since they are vacuum properties and therefore 
universal, they can be extracted from phenomenological analyses of hadron masses, as it is 



customary done in QCD sum rules 



The condensates are expected to vanish in the limit 



of very high temperature or chemical potential. However, it has been suggested that they 
may survive after the deconfinement transition both in the high temperature 



llO| and 

in the high chemical potential cases For our purposes the relevant gluon condensates 
are those of dimension four 0], (0|^F^|0) (= (F^)), and of dimension two Q, (ol^f^A^IO) 

We shall derive an equation of state (EOS) for the cold QGP, which may be useful for 
calculations of stellar structure. Our EOS can be considered an improved version of the EOS 
of the MIT bag model, which contains both the non-perturbative effects coming from the 
residual gluon condensates and the perturbative effects coming from the hard gluons, which 
are enhanced by the high quark density. As it will be seen, the effect of the condensates is 
to soften the EOS whereas the hard gluons significantly increase the energy density and the 
pressure. 



II. THE EQUATION OF STATE 



In this section we develop a mean field approximation for QCD, extending previous works 
along the same line 13l-ll7l|. The Lagrangian density of QCD is given by: 



QCD 



4 



q=l 



^1 



(1) 



where 

The summation on q runs over all quark flavors, is the mass of the quark of flavor g, 
i and j are the color indices of the quarks, are the SU(3) generators and /"'"^ are the 
SU(3) antisymmetric structure constants. For simplicity we will consider only light quarks 
with the same mass m. Moreover, we will drop the summation and consider only one flavor. 



At the end of our calculation the number of flavors will be recovered. Following 13|, llJ], we 
shall start writting the gluon field as: 



where A"''^ and a"'^ are the low ("soft") and high ("hard") momentum components of the 
gluon field respectively. The former will be responsible for the long range and low momentum 
transfer, non-perturbative processes whereas the latter will be relevant in the short distance 
perturbative processes. The field decomposition made above requires the choice of an energy 
scale defining the frontier between soft and hard. This energy scale, E, lies in the range 
^QCD < E < 1 GeV. In principle, the dependence of the results on this choice can be 
studied with the renormalization group techniques. Accurate results would also require the 
knowledge of the scale dependence of the in-medium gluon condensates, which in our case 
is poor. Therefore, in order to keep the simplicity of our approach, we will not specify the 
separation scale and will assume that A"-^ represents the soft modes which populate the 
vacuum and a""^ represents the modes for which the running coupling constant is small. 
Inserting ([3]) into ([2]) we obtain: 

+ gr^^A^^'a"'' + gf^'^a^^'A^'' (4) 

In the above expression the coupling is running and is large (small) when attached to A"-^ 
[a""^). The mixed terms, such as gf°''"^a'''^A'^'^, are assumed to be dominated by the large 
couplings. 



A. The mean field approximation 



In a cold quark gluon plasma the density is much larger than the ordinary nuclear matter 
density. These high densities imply a very large number of sources of the gluon field. With 
intense sources the bosonic fields tend to have large occupation numbers at all energy levels, 
and therefore they can be treated as classical fields. This is the famous approximation for 



bosonic fields used in relativistic mean field models of nuclear matter [18|. It has been 
applied to QCD in the past |l5| and amounts to assume that the "hard" gluon field, a^, is 
simply a function of the coordinates 18| : 

< = "o V (5) 
n fact, for cold nuclear matter, it is further assumed that aj^ is constant in space and time 

d,a; = (6) 



As a consequence of this approximation, the term gf^'^'^a^aQ will vanishe because of the color 
symmetry. We also assume that the soft gluon field A"-^ is independent of position and time 
and thus: 

d^A""" = (7) 

Substituting ([5]), ([6]) and d?]) into (gD we have F"^^ = ^/'^^^(A^'^A^^ + A^/^aJ^^^" + a^^^'M^^). 
Inserting this into ([T]), the QCD Lagrangian simplifies to: 

g2 Jabc Jade 
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+AlAlA^>'al5'^'' + Aj.Al^a^^O'^A"'^ + Alal5Q,A'^^' A'" + alS^^AlA'^'' A'" 



(8) 



We shall now replace the soft gluon field A^^ and its powers by the corresponding expectation 
values in the cold QGP. The product of four fields in the first line of the above equation can 

nn 

be related to the gluon condensate through the relations pj|, 1141 ]: 



{A^AlA^P A^"^) 



(32) (34) 



m irab jrcd , p ri irac roa , ri p sraa ri 



and 



_ / pa^lu p \ 



where the constant b is given by: 



g^^TT 
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4(34)' 



jpafiu pa ^ 



(9) 



(10) 



(11) 



In the second and fourth lines of ([8]) we have odd powers of A""^ which have vanishing 
expectation values: 

{A^^'A^'^A'P) = (12) 

(A"^) = (13) 

we have the hard gluon mass terms. The expectation value of two 



In the third line of 
soft fields reads 
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{A^PA^") 



'^2~^° 



(14) 



The {g'^A^) gluon condensate alone is not gauge invariant. While this might be a problem in 
other contexts, here it is not because it appears always multiplied by other powers of gluon 



fields, forming gauge invariant objects. The {g^A^) condensate is associated [is], [l^ with a 
dynamical gluon mass: 

In spite of the recent progress in the field, still very little is known about {A^) at finite 



(and high) density. In our approach, as in [l9|, we have {g^A^) < in (fT4|) so is positive. 



Using expressions ([9]), ( fT2l) . ( IT3l) . ( fT4|) and ( lT5l) in ([H]) we arrive at the following effective 
Lagrangian: 
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B. Pressure and energy density 

From the Lagrangian (fT6|) we can derive the equations of motion: 

ma'al = -gp'' (17) 

(t S,j Yd, + 9l'^< - ^) =0 (18) 
where p° is the temporal component of the color vector current given by: 

= Al'T^.i^j (19) 

From the Lagrangain we can obtain the energy-momentum tensor and the energy density of 
the system through: 

e =< Too > (20) 
In the present case the energy-momentum tensor is given simply by: 

and consequently: 

d£ 

which, with the use of f lT6|) gives: 

2 

e = ^iJ^J%^oi^^) - goo{ - b4>o' + ^ao"«o + + gi'^a^ - (5,,m)^]} (23) 

Using f lTSj) in the expression above we find 

2 

e = b(Po' - ^«o«S + ^^^7'{^o^^) (24) 



Multiplying (fT8|) by ipi from the left we find: 

i^idoi'i) = ijli-id ■ V + 7°m)^^, - (25) 



From the usual Dirac theory applied to the study of nuclear matter we have [18| : 

ipli-ia ■ V + 7°m)^i = 3^ [ ^ dk k'^Jp + m^ (26) 
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In the last two expressions we have: 

'=(;: 

where a are the standard 2x2 Pauli matrices, the unit entries in 7° are 2x2 unit matrices 
and 7q is the quark degeneracy factor 7q = 2(spin) x 3(fiavor). The sum over all the color 
states was already performed and resulted in the pre-factor 3 in the expression above, kp is 
the Fermi momentum defined by the quark number density p: 



which gives: 



P^^M (27) 



In the above expression |A^) denotes a state with N quarks. Inserting (126|) into (l25l) and 
then ([25]) into ^ we find: 

'-a^a^o - gp^al + 3^ f ^ dkk'^ \lp + m'^ (28) 
2vr Jq 



Using ( !T7|) we can eliminate the field in the above expression: 



We can relate the color c harg e density and the quark number density p. To do this we 



shall use the notation of 20|] and write the quark spinor as ipi = ipCi, where q is a color 
vector. We have: 

pV^ = {^Pa'mj.m^'nM = {4^^3MnA) = (cj7;^.c,)^V(47^.V/)^V = 3p^ 

(30) 



where we used the relations ^^^p = p and {(^^T^-Cj){(^^T^iCi) = 3. Performing the momentum 
integral we arrive at the final expression for the energy density: 



> + 600 + ^TT^i + a 



2771^2 y ^ 27r2 1 4 

4 



m 
— In 



I "1 TTl^ 1 

kF+ykp^ + TO^ H /nfm^) > 

J 16 J 



(31) 



The pressure is given by 

P = ^ < T,, > (32) 
Repeating the same steps mentioned before and using: 

VW, = 3-^ / A.|^j^}=3^ r<i«4-ji=} (33) 

we arrive at: 

P = ^«o«o - Ho' + ^ Tdkel I (34) 

Performing the momentum integral, using (fT7|) and the relation for p*^ and the quark number 
density p in fl5^ we obtain the final expression for the pressure: 



3^2 \ ( kp^Vkp^ + ?,m^kF\/kF^ + 

: ) P - 000 + : 



Wj^ -^'^^ ' 2.rM 4 8 



3m^ , 
H /ra 



fc^ + + m2 1 - ^/r2(m2) \ (35) 

J 16 J 



The speed of sound Cg is given by: 



2 {oa\ 
cs =g-^ (36) 



In the expressions above, g is small, since it comes always from the coupling between the 
hard gluons and the quarks. The large coupling is contained in the constants b and mo- 

Both fl3T]) and fl35|) have three terms. The first term, proportional to p^, comes from 
the purely hard gluonic term appearing in the Lagrangian and from the hard gluon term 
appearing in the quark equation of motion. The second term, proportional to comes 
exclusively from the soft gluon terms and it has opposite signs in the energy and in the 
pressure. This is precisely the behavior of the bag constant term in the MIT bag model 
which has the same origin. The third term comes from the quarks. In short, we can say 
the both the energy density and the pressure are the sum of three contributions: the hard 
gluons, the soft gluons and the quarks. 
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III. NUMERICAL RESULTS AND DISCUSSION 



We now compare our results (!3T]) . (!35|) and (!36l) . with the corresponding results obtained 



with the MIT bag model for a gas of quarks at zero temperature 



eo={^W/'pB'/' + B (37) 



and 

Po = ^(^)^^/W/^-S (38) 

and 

CO = ^ = 3 (39) 
We choose B = 110 MeV fm"^, which lies in the range (50 < B < 200 MeV fm~^) used in 



calculations of stellar structure 



22 



2|-|24|. For the comparison we must rewrite fl271) . (|3T1) and 
(l35l) as functions of the baryon density, which is ps = |p- 

If we neglect the gluonic terms and choose the quark mass m to be zero in ( 1311) and ( 135|) 
we can show that they coincide with ( 1371) and (1381) with B = 0. In this limit, our model 
reduces to the MIT bag model. 

We next consider the MIT bag model with finite B and our model with massless quarks 
and soft gluons but no hard gluons. This comparison is meaningful because with these 
ingredients both models describe the dynamics of free quarks under the influence of a soft 
gluon background. In this case we can identify our gluonic term with the gluonic component 
of the MIT bag model, represented by the bag constant. We then obtain an expression for 
the bag constant in terms of the gluon condensate: 

Bqcd = b4 = {\f^^''FP (40) 

where in the last equality we have used J 10]) and flTTj) . The above relation has been found in 
previous works, such as, for example, |9|. Fixing B and choosing a reasonable value of the 
coupling of the soft gluons, appearing in ( ITOj) we can infer the value of the dimension four 
condensate, < F"^ >, in the deconfined phase. For B = 110 MeV fm~^ and g = 2.7 (which 
would correspond to as = 5'^/47r = 0.6) we find: 

(F^) = (^F'/^-F" ) = 0.0006 GeV^ (41) 

71 
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In the lack of knowledge of the in-medium dimension two condensate, we use the factor- 



ization hypothesis, which, in the notation of Refs. 



13| and jl^, implies the choice /io 



As a consequence, (jH]), (HI?]) . (ITll) and (HM are related and we obtain: 



{g'A^) = -y^llKM![!(F2) = _o.3 GeV^ (42) 
which corresponds to a dynamical mass of rriG = 290 GeV. This number is consistent with 



the values quoted in recent works [25|-|27|, which lie in the range 200 < rriG < 600 MeV. 
Finally, the numerical evaluation of f l3T|) . fl35|) and f p6|l requires the choice of g, the coupling 
of the hard gluons, and of the quark mass, m. We choose them tohe g = 0.35 (corresponding 
to as = g^/^TT = 0.01) and m = 0.02 GeV. 

In Fig. 1 we show the energy density, pressure and speed of sound obtained with (13T1) . 
(153]) and fl5B]) divided by the corresponding MIT values: Eq, po and Cq. We observe that, 
for this set of parameters, our EOS is harder than the MIT one. This can also be seen in 
the plot of the pressure as a function of the energy density, shown in Fig. 2. In the same 
range of baryon densities, we have more energy, much more pressure and consequently a 
larger speed of sound. This behavior can be attributed to the first term of the equations 
( pTl) and ( 135|) . which comes from the hard gluons. This term is exactly the same both in 
(13T1) and (l35l) and in the limit of high densities becomes dominant yielding p ~ e and hence 

-> 1. Physically, this term represents the perturbative corrections to the MIT approach. 
Since the quark density is extremely large, even in the weak coupling regime (typical of the 
hard gluons) the field is intense. A similar situation occurs in the color glass condensate 
(CGC). In that context, a proton (or nucleus) boosted to very high energies becomes the 
source of intense gluon fields generated in the weak coupling regime. Also in that case 
semi-classical methods were applied to study this gluonic system. 

In Fig. 3 we plot the energy density fl3T]) (upper panel) and the pressure fl35|) (lower 
panel) as a function of the baryon density pb- We take pb always starting at 1.5 fm~'^. 
We can observe that the quarks and hard gluons give the dominant contributions both to 
the energy and to the pressure. Looking at the pressure we see that the hard gluons give a 
repulsive contribution whereas the soft gluon contribution is attractive. It is interesting to 



see that our curves follow very closely those of Refs. 23|] and 2^], computed with slightly 
different versions of the MIT bag model. 

In Fig. 4 we show the EOS for different choices of the condensates, which are now treated 
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as independent from each other. In the upper panel we fix (F^) and vary {g^A^), starting 
from the central value —0.3 GeV^ and increasing its magnitude. In the lower panel we 
perform the complementar study keeping (g^A^) and increasing the magnitude of {F"^). As 
it can be seen, increasing the condensates reduces the pressure and, in the case of {g'^A'^), 
softens the equation of state. This behavior could be anticipated from Eqs. fl3ip . (1351) and 
from equation of motion (fT71) . Indeed, keeping fixed the coupling and the quark density, 
when we increase the gluon mass, the field becomes weaker. In a more accurate treatment, 
with the inclusion of spatial inhomogeneities, the equation of motion (fT7j) would contain a 
Laplacian term and its solution would show a Yukawa behavior, with the mass mc controlling 
the screening of the field Oq. 

In Fig. 5 we show the energy per particle as a function of the baryon density for different 
values of the gluon condensates. As in the previous figure, in the upper panel we fix (F^) 
and vary {g^A"^). Increasing {g'^A^) the energy per particle grows slower with baryon density. 
The system becomes more compressible. In the lower panel we keep {g'^A^) fixed and increase 
the magnitude of {F"^). Increasing (F^) leads, as before, to a more compressible system but 
the total energy is now larger. For the central values of (F^) and {g'^A^) we obtain values of 



e/ pb which are compatible with those found in Ref. [28j for equivalent baryon densities. As 
it can be seen in all curves, the energy per particle is always much larger than the nucleon 
mass (939 MeV) and hence the system under consideration can decay into nuclear matter. 

To summarize, we have derived an equation of state for the cold QGP, which may be useful 
for calculations of stellar structure. The derivation is simple and based on three assumptions: 
i) decomposition of the gluon field into soft and hard components; ii) replacement of the 
soft gluon fields by their expectation values ( "in-medium condensates" ) and iii) replacement 
of the hard gluon fields by their mean-field (classical) values. Our EOS can be considered 
an improved version of the EOS of the MIT bag model, which contains both the non- 
perturbative effects coming from the residual gluon condensates and the perturbative effects 
coming from the hard gluons, which are enhanced by the high quark density. It is reassuring 
to observe that our EOS has the correct limits, where we recover the MIT bag model results. 
The parameters are the usual ones in QCD calculations: couplings, masses and condensates. 
The effect of the condensates is to soften the EOS whereas the hard gluons significantly 
increase the energy density and the pressure. 
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FIG. 1: Energy density, pressure and speed of sound, as functions of the baryon density, divided 
by the corresponding MIT values: Eq, po and cq. 
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